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Abstract 

In this paper we study the evolution operator of a time-dependent Hamiltonian in 
the three level system. The evolution operator is based on SU(3) and its dimension 
is 8, so we obtain three complex Riccati differential equations interacting with one 
another (which have been obtained by Fujii and Oike) and two real phase equations. 
This is a canonical form of the evolution operator. 

1 Introduction 

In this paper we treat a finite dimensional quantum model and study the evolution operator 
of a time-dependent Hamiltonian from a geometrical point of view. Then we must solve the 
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time-dependent Schrodinger equation, which is a very hard task. Our method is based on 
one in PQ, and 0. 

Let us start by setting the problem. By Ho(n; C) we show the set of all n x n-hermitian 
matrices with zero trace. The general form of H is given by 
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with h^t) + h 2 it) + h 3 (t) + ■■■ + h^it) + Kit) = 0. 



Under this Hamiltonian we want to find the evolution operator U = U(t) satisfying the 
Schrodinger equation 

iU = i^£ = HU (2) 

dt y ' 

where we have set % = 1 for simplicity. The wave function \I/ = \l/(t) is of course related to 
the evolution operator like ^Sf(t) = Z7(i)^(0), i7(0) — E (: identity). However, it is almost 
impossible to solve (j2J) exactly except for a few examples. 
We decompose U(t) into two unitary parts 



U(t) = U^Z{t))U 2 {t) 



(3) 



where Z is a local coordinate of some symmetric space like Grassmann manifolds or (gener- 
alized) flag manifolds, see for example jl], [5], [6] and [7]. In the case of flag manifolds this 
decomposition is based on the approximation 



SU(n) « {^(n)/^!) 71 - 1 } x L^l)*" 1 . 



(4) 



In the following sections we treat the Bloch sphere SU(2)/U{1) = CP 1 = S 2 and the 
flag manifold S'f/(3)/[/(l) 2 as interesting examples. 



U2(t) is a small unitary part and in the paper it is set as the phase part 



oi<h(t) 



\ 



i<j> n (t) 



; 0i(t) + 2 (t) + --- + 0„(t) = 0. (5) 
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Then the equation (j2J) is reduced to a complicated combination of equations on Z(t) = 
(zij(t)) and {4>k(t)}. We want to call it a canonical form in our sense. Unfortunately, it is 
very hard to write down the general case explicitly. 

In this paper we transform the equation §2§ into a canonical form in the case of n = 2 and 
n = 3. Namely, since dim^ SU(2)=3 for n = 2 we obtain one complex Riccati differential 
equation and one real phase equation. See for example [1] and its references. 

On the other hand, since dim^ S77(3)=8 for n = 3 we obtain three complex Riccati 
differential equations interacting with one another (which have been obtained by Fujii and 
Oike [7]) and two real phase equations. 



2 Two Level System 

In this section we treat the two level system (n = 2) in detail. The Hamiltonian is 

, h(t) v(t) , 

(6) 



from ([T]). From the fact 



v(t) -h(t) 



SU(2) « x 1/(1) <* S 2 x 17(1) 



(see [7J), U(t) can be parametrized as 

1 / 1 -z(t) \ ( e**® 
U{t) = Ui{Z{t))U*{t) = | (7) 

V l + \ z (t)\ 2 \ z(t) 1 J \ e-^W 

From the equation (j2J) we have only to calculate 

H = iUU ] =iUU~ 1 . 
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Some calculation gives 

ov* = _#±ii l2+ ^u(° -*)( 1 M+^f 1 
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and from (181) we have 
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or 
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From this it is not difficult to show 



j^{z + 2i<fiz) = v, 

i# - m - n 2 )} = -/» 

z + 2i0z = — if (1 + |z| 2 ), I — 2i(f)z = iv(l + |z| 2 ), 
zz ~ z% - i<p{l - \z\ 2 ) = ih(l + |z| 2 ). 



z + 2i<fiz = —iv(l + |z| 2 ), 

i v z+vz+2h 



and finally we obtain 



z - i(vz 2 + 2hz - v) = 0, (9) 
0=- ^ + ^ + 2/ \ (10) 

The equation (jUJ) is a Riccati equation and (II Op is a phase equation. This is our canonical 
form of (jHJ). We believe that our derivation is smarter than that of pQ. 

A comment is in order. Unfortunately, it is impossible to solve fl9]) exactly, so we need 
some numerical computational method. 



3 Three Level System 

In this section we treat the three level system {n — 3) in detail, which is the aim of the 
paper. However, the calculation is not easy compared to the two level one. 



The Hamiltonian is 
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where 
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and 



and 



Al = 1 + | x (t)| 2 + \y{t)\\ A 2 = 1 + \z(t)\ 2 + \x(t)z(t) - y(t)\ 2 (14) 
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This form is convenient in the latter calculation, though it is a variant of that of [7], [8], [9]. 
Moreover, for the latter convenience we set 
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Then U(t) in (|T2j) can be written as 

U(t) = V(t)D A (t)U 2 (t). 

In the following we omit t from all variables {x(t) — > x, etc) for simplicity. From the 
equation (T5]) we have 

i (VD A U 2 + VD A U 2 + VD A U 2 \ = HVD A U 2 

and 

i (V + VD A D^ + VU 2 U^ = HV (18) 
because D A and U 2 are diagonal. This form is better as shown in the following calculation. 



A comment is in order. From flTBl we have a clear form 

h = i (vv- 1 + vb^-^v- 1 + vii 2 u 2 - V" 1 

However, we don't recommend readers to calculate the right hand side because the calculation 
becomes very complicated. 

The calculation is truly complicated. The (11), (21) and (31)-components of the matrix 
equation ffT8l read 

.. «Aj • 

ix — x — (pix = Vi + h 2 x + v 3 y, 

2 Zai 

iAi • 

0i = h 1 +v 1 x + v 2 y. 

2 Zai 

From these equations we have 

ix = (v i + h 2 x + v 3 y) - (hi + v x x + v 2 y)x 

= v\ + (h 2 - hi)x - vix 2 + - v 2 xy (19) 

and 

«y = (v 2 + v 3 x + h 3 y) - (hi + v x x + v 2 y)y 

= v 2 + (h 3 - hi)y - v 2 y 2 + v 3 x - v x xy (20) 
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and 

X lzAl 
Next, we must calculate the term 



(h 1 +v 1 x + v 2 y). (21) 



i/S.i = i— (1 + \x\ 2 + |y| 2 ) = i(xx + xx + yy + yy). 



dt 

Then from (USD and fl20l) it is not difficult to show 



iAi = {(/ii + vix + v 2 y) - (/ii + vix + v 2 y)}Ai, 
so that we finally obtain 

• (hi + v\x + v 2 y) + (hi + vix + v 2 y) 



(22) 



from (12T|). 

The (12), (22) and (32)-components of the matrix equation fll8p read 
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From these equations it is not difficult to show 



iz = v 3 + h 3 z - (h 2 + v 3 z)z + (xz - y){v l + v 2 z) 

= v 3 + (h 3 - h 2 )z - v 3 z 2 + (xz - y)(v! + v 2 z) (23) 
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(/i2 + ^3-2) + ^(Vi + V 2 z). 



iAi 

— = (h x + vix + v 2 y) - {hi + fix + v 2 y) 

the remaining one is to calculate A2/A2. However, it is very troublesome. 
From 

iA 2 = i— (1 + \z\ 2 + \xz-y\ 2 ) 
at 



i I zz + zz + -f- (xz — y) (xz — y) + (xz — y)-r (xz — y) \ 
I at at J 



long calculation gives 

iA 2 



-t> 3 z + f 3 z + u 2 (x2; — y) — v 2 (xz - y) 



A 2 

by use of ffT9]) . f )20|) and f l23|) . Therefore we finally obtain 

• (— h 2 - v 3 z + ViX + v 2 xz) + {—h 2 - v 3 z + vix + v 2 xz) 

02 = 

from (1241). 

Let us summarize the result. Our canonical form is 

ix = vi + (h 2 — hi)x — v\x 2 + v 3 y — v 2 xy, 
iy = v 2 + (h 3 - h{)y - v 2 y 2 + v 3 x-v x xy, 
iz = v 3 + (h 3 - h 2 )z - v 3 z 2 + {xz - y)(vi + v 2 z) 

■ (hx + ViX + v 2 y) + {h\ + Vix + v 2 y) 

0i = g ' 

• (—h 2 - v 3 z + v x x + v 2 xz) + (—h 2 - v 3 z + Vix + v 2 xz) 



and 



where hi + h 2 + h 3 = 0. 

A comment is in order. We have a set of complex Riccati differential equations interacting 
with one another ([7j) and real phase equations. Unfortunately, to solve them exactly is 
impossible, so we must use some numerical computational method. 



4 Discussion 

In this paper we obtained a canonical form of the evolution operator of a time-dependent 
Hamiltonian in the three level system (n = 3), namely three complex Riccati differential 
equations interacting with one another and two real phase equations. For the case of n > 4 
it is very hard to calculate, so we will leave calculation to (young) readers. 

The result looks good and may be applied to a wide class of problems in Quantum 
PhysicJi] or Mathematical Physics. Further work will be needed. 
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